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We investigate the capability of directional detectors to probe neutralino dark matter in the Min- 
imal Supersymmetric Standard Model and the Next-to-Minimal Supersymmetric Standard Model 
with parameters defined at the weak scale. We show that directional detectors such as the future 
MIMAC detector will probe spin dependent dark matter scattering on nucleons that are beyond 
the reach of current spin independent detectors. The complementarity between indirect searches, in 
particular using gamma rays from dwarf spheroidal galaxies, spin dependent and spin independent 
direct search techniques is emphasized. We comment on the impact of the negative results on squark 
searches at the LHC. Finally, we investigate how the fundamental parameters of the models can be 
constrained in the event of a dark matter signal. 



PACS numbers: 95.35. +d, 14.80.Ly 

I. INTRODUCTION 

The existence of cold Dark Matter (DM) is inferred 
from a large number of astrophysical observations at var- 
ious scales. Locally, the large discrepancy between New- 
ton's law of gravitation and the observed rotation curves 
indicates that spiral galaxies should be embedded in a 
DM halo [U[2]. 

Candidates for this class of to-be-discovered particles nat- 
urally arise from extensions of the standard model of par- 
ticle physics (e.g Supersymmetry), as long as an ad-hoc 
discrete symmetry is invoked to preserve the proton sta- 
bility. Indeed, the lightest particle of the secluded sector 
(in Supersymmetry: the lightest supersymmetric particle 
(LSP)) becomes stable and can be a good candidate for 
the cold DM present in the Universe especially if it is a 
weakly interacting massive particle (WIMP). Numerous 
DM candidates have been proposed in extensions of the 
standard model and include fermions, scalars or gauge 
bosons [3]. 

Tremendous experimental efforts on a host of techniques 
have been made in the field of WIMP detection, the main 
experimental issue being the WIMP-background discrim- 
ination. Directional detection of galactic DM has been 
proposed as a powerful tool to identify genuine WIMP 
events as such [4]. Recent studies have shown that, 
within the framework of dedicated statistical data analy- 
sis, a low exposure directional detector could lead either 
to a high significance discovery of galactic DM [5H7] or 
to a conclusive exclusion [8 , depending on the value of 
the unknown WIMP-nucleon cross section. In the case of 
a high significance detection, it would also be possible to 
go further and to constrain the WIMP properties, both 
from particle physics (mass, cross section) and galactic 
halo (velocity dispersions), within the framework of a 
model independent analysis [9]. To achieve this would 



require a rather high spin dependent (SD) cross section 
(of the order of l(T 4 pb) and a low WIMP mass, com- 
parable to the target nucleus mass, as a matter of fact 
below the electroweak scale. Most projects of directional 
detection [10] are low pressure Time Projection Cham- 
bers (50-100 mbar). It follows that there is an obvious 
size limitation, which in turn implies that directional de- 
tectors cannot scale up to ton-scale experiment as most 
of the current direct detection projects [TTH23] . Then, 
to be competitive and complementary with planned and 
existing direct detectors, directional detection should fo- 
cus on SD interaction for which a large fraction of the 
parameter space could be probed by planned directional 
experiments (see SecjlT]). 

DM models which can best be probed by directional SD 
detection are the ones with Majorana particles such as 
the neutralino in the Minimal Supersymmetric Standard 
Model (MSSM) or its extensions. Not only can the neu- 
tralino be light (below 100 GeV) but also its Majorana 
nature implies that different processes contribute to the 
SD or spin independent (SI) interaction, namely the Z 
exchange contributes only to the SD cross section. Both 
types of detection modes are therefore complementary. 
Vector DM candidates can also lead to a signal in the SD 
mode which is not directly correlated with the SI mode. 
However in most models studied previously (for exam- 
ple UED or little Higgs [24] [25]) the DM tends to be 
heavier than 100 GeV (even heavier than 1 TeV in the 
minimal UED case [26]) making it difficult to extract 
the DM properties for direct detection. Models with 
Dirac fermion DM get a contribution from Z exchange 
for both SI and SD, these DM candidates are therefore 
best probed by SI where one can take advantage of the 
coherence effect. 

We will therefore consider only models with a Majo- 
rana DM particle, more specifically the neutralino in the 
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MSSM and in one of its extensions, the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM) [27]. The 
parameter space of these models consistent with collider 
physics, precision measurements, DM relic density as well 
as DM direct searches were examined in several studies 
assuming either some GUT scale relations among the fun- 
damental parameters [28-36] or the more general case of 
parameters defined at the electroweak scale [37U42] . The 
potential of indirect detection for probing these models 
was examined for example in [43 -49 and the complemen- 
tarity between detection techniques were investigated 
in [501153] . In particular, in a previous study [54] [55] 
the allowed parameter space of the MSSM and of the 
NMSSM with light neutralinos (up to 50 GeV) was ex- 
plored using a MCMC approach. Furthermore indirect 
detection of DM through the flux of photons from dwarf 
spheroidal galaxies, as well as direct detection limits from 
XENON100 were shown to be complementary and to fur- 
ther constrain the parameter space. 
In this study we first determine the parameter space of 
the MSSM and NMSSM that can be explored by direc- 
tional SD interactions with the future MIMAC detec- 
tor [56]. For this we consider in both models the whole 
spectrum of masses for neutralinos although we use a re- 
stricted set of free parameters at the electroweak scale. 
We then compare this with limits from SI and indirect 
detection searches, namely gamma-rays. Finally, we in- 
vestigate how the fundamental parameters of the models 
can be constrained in the event of a DM signal. 

The paper is organized as follows. In section [II] 
we discuss the potential of SD directional detector in 



probing DM. In section III we discuss the expectation 
for the SD cross section. The parameter space of the 
model considered and the scanning method used are 
presented in Section IV Our results for the SI/SD direct 
searches as well as for the flux of photons relevant for 
indirect searches are presented in Section |V| highlighting 
the complementarity between different techniques. The 
impact of directional detectors on probing the parameter 
space of the model is also addressed while the impact 
of the LHC on these results is sketched in section IVII 
An investigation on how the observation of a signal 
could be used to determine some parameters of the 
supersymmetric model is carried in section \VT\\ 



spectrum d 2 R/dE r dVt r down to the energy threshold. It 
is worth emphasizing that it is the lowest energy at which 
both the track and the energy can be retrieved which is 
the key experimental issue for directional detection. It 
follows that, to maximize the track length, the pressure 
of the gaseous detector must be as low as possible, lead- 
ing to rather small detector masses as the volume cannot 
be arbitrarly large. One may then come to the conclu- 
sion that the directional detection strategy should focus 
on SD interactions to be competitive with planned and 
existing direct detectors [TTH23] . 

Then, the ideal directional target is a nucleus with non- 
vanishing spin. Leading candidates include: 1 H, 3 He and 
19 F which has been early suggested as a golden target 
for SD DM searches [62 . CF4 is indeed planned as a 
sensitive medium for most upcoming directional detec- 
tors [10]. 

In the following, as a working example, we present the 
case for a low exposure (30 kg. year) CF4 TPC, operated 
at low pressure and allowing 3D track reconstruction, 
with sense recognition down to 5 keV. Such performance 
is taken as the ultimate limit for a directional detector. A 
complete overview of the effect of the main experimen- 
tal issues of directional detectors, such as background 
contamination, energy threshold, sense recognition effi- 
ciency, angular and energy resolution on exclusion limits 
and discovery potential, are presented in [6j[8]. 



B. Directional detection 

1. Directional event rate and astrophysical inputs 

Directional detection depends crucially on the local 
WIMP velocity distribution [63-65 . The isothermal 
sphere halo model is often considered but it is worth go- 
ing beyond this standard paradigm when trying to ac- 
count for all astrophysical uncertainties. The multivari- 
ate Gaussian WIMP velocity distribution corresponds to 
the generalization of the standard isothermal sphere with 
a density profile p(r) oc 1/r 2 , leading to a smooth WIMP 
velocity distribution, a flat rotation curve and no sub- 
structure. The WIMP velocity distribution in the labo- 
ratory frame is given by, 



II. DIRECTIONAL DETECTION FRAMEWORK 

A. Directional detectors 

There is a worldwide effort toward the development of 
a large TPC (Time Projection Chamber) devoted to di- 
rectional detection [10] and all current projects [56H6T] 
face common experimental challenges and share a unique 
goal: the simultaneous measurement of the energy (E r ) 
and the direction of the 3D track (Q r ) of low energy 
recoils, thus allowing to evaluate the double-differential 



m = 



1 



Btt 3 det<r2)i/2 



exp 



(v-Vq) a v (v-Vq) 



(1) 

where cr v = diag^, <j y , a z ] is the velocity dispersion ten- 
sor assumed to be diagonal in the Galactic rest frame (£, 
z) and v Q is the Sun's velocity vector with respect to 
the Galactic rest frame. When neglecting the Sun pecu- 
liar velocity and the Earth orbital velocity about the Sun, 
Vq corresponds to the detector velocity in the Galactic 
rest frame and is taken to be Vq = 220 km.s -1 along 
the y axis pointing toward the constellation Cygnus at 
(£q = 90°, b@ = 0°), where £ and b are the Galac- 
tic latitude and longitude. This way, we can consider 
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the three velocity dispersions along the three axis as nui- 
sance parameters in order to take into account the effect 
of anisotropic WIMP velocity distribution when deriving 
the discovery potential of upcoming directional detectors. 
The anisotropy is defined by the f3 parameter as: 



/3 = 1 



2a? 



(2) 



Hence, j3 = corresponds to an isotropic WIMP velocity 
distribution while /3 < corresponds to a tangential 
anisotropy and f3 > to a radial anisotropy. Using 
the parametrization of the three velocity dispersions 
used in this study (see table |l|, we are considering 
WIMP velocity distributions with /3 = ± 0.25 (68% 
C.L.). This range is compatible with recent results from 
N-Body simulations which found j3 = — 0.4 [73] and 
observations which found a non vanishing j3 parameter 
in the solar neighbourhood [74] [75] . 



The directional recoil rate is given by [66 
d 2 R 



dE r dQ r Ai:m x m* 
with rn x the WIMP mass. 



F z (E r )f(v min ,q), 



(3) 



m r the WIMP-nucleus 
reduced mass, po the local DM density, ctq the WIMP- 
nucleus elastic scattering cross section, F(E r ) the form 
factor (using the axial expression from [67]), v m i n the 
minimal WIMP velocity required to produce a nuclear 
recoil of energy E r and q the direction of the recoil 
momentum. Finally, f(v m i n ,q) is the three-dimensional 
Radon transform of the WIMP velocity distribution 
/(#), see [66 for more details. 

As one can see from eq. [3] the directional rate is directly 
proportional to the local DM density at Solar radius 
(po) which is also subject to important uncertainties. 
Following [6 we consider po as a nuisance parameter, 
using po = 0.3 ± 0.1 GeV/cm 3 . We used a mean value 
of 0.3 GeV/cm 3 for the sake of comparison between the 
different direct searches experiments. 
The last astrophysical uncertainties to be considered 
when deriving the discovery potential of upcoming di- 
rectional detectors is the the velocity of the Solar system 
orbit's in the Galactic rest frame taken as Vq = 220 ± 30 
km/s. 

It is worth noticing that other astrophysical uncer- 
tainties like the escape velocity could be taken into 
account when assessing the sensitivity of a given direct 
detection experiment. However, the escape velocity is 
not considered in this study as we are dealing with a low 
mass target material ( 19 F) associated to a low energy 
threshold. Indeed, as the minimal speed required to 
produce a 5 keV recoil energy is about v m i n = 130 km/s 
for a WIMP mass of 100 GeV/c 2 , while the mean WIMP 
velocity in an earth based detector is about 300 km/s, 
we can deduce that the effect of a finite escape velocity 
will be negligible. 



Nuisance parameters 


Gaussian parametrization 


nn rripV/r 2 /rm 3 l 




VQ [KII1/SJ 


OOf) _|_ Of) 


<j x [km/s] 


220/^2 ±20 


g v [km/s] 


220/^ ±20 


g z [km/s] 


220/^ ±20 



TABLE I: Gaussian parametrization (mean and standard de- 
viation) of the different astrophysical nuisance parameters. 



As a conclusion, accounting for uncertainties on the 
astrophysical parameters is a step beyond the "standard 
DM model" , i.e. isotropic isothermal DM halo, with fixed 
values of density and Sun's circular velocity. Evaluating 
the properties of the DM halo is indeed still a subject of 
debates. The numerical values and uncertainties corre- 
sponding to these astrophysical inputs are given in table [I] 
(see [6] for a detailed discussion). 
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FIG. 1: Spin dependent cross section on proton (cm 2 ) as a 
function of the WIMP mass (GeV) in the case of the pure- 
proton approximation [78] showing the sensitivity of a forth- 
coming 30 kg. year directional experiment (solid red line). It 
is denned as the minimal cross-section required to obtain a 
DM discovery with a significance greater than the 3a level at 
90% C.L. [6]. For the same exposure, the dotted black line 
presents the projected exlusion limit [8 in the background 
free case. The red star is the input value of the benchmark 
model and the black contour is the 68% contour level ob- 
tained with the MCMC analysis [9]. The exclusion limit from 
COUPP-2010 [11] (black dotted line) is also presented. 
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2. Dark Matter properties inferred from directional 
detection 

As highlighted by several recent studies [5j [7H9] , direc- 
tional detection may be a powerful tool to discriminate 
between the DM signal and the background one. Indeed, 
the correlation between the main incoming direction of 
the recoiling events in galactic coordinates with the So- 
lar system orbit's around the Galactic Center has been 
shown to be a strong and convincing proof in favor of a 
DM detection. This kind of discrimination is even more 
relevant when considering the recent results from DAMA, 
CoGeNT and CRESST experiments which have observed 
candidate events which origin are difficult to assess with- 
out directional information. 

At first, one may think of using directional detection to 
set exclusion limits. Several methods have been pro- 
posed [8j [79]. Figure [I] presents the expected rejection 
limit for a 30 kg. year CF4 directional detector, in the 
background-free case [8 for which a standard limit has 
been derived using the classical Poisson statistics associ- 
ated to a observed event. It allows to reach ~ 10 -6 pb, 
noticing that for highly background-contaminated data 
(~ 10 kg -1 . year -1 ) the result would be about an order 
of magnitude higher. 

However, directional detection may be used to go be- 
yond the standard exclusion limit strategy. Indeed, it 
may allow to discover DM [5j [7], the proof of discov- 
ery being the fact that the signal points to the direction 
of the constellation Cygnus (to which the solar system's 
velocity vector is pointing). Hence, the goal is to iden- 
tify a genuine WIMP signal as such. Using a frequentist 
profile likelihood ratio test statistics, and taking into ac- 
count astrophysical and experimental uncertainties, one 
can determine the sensitivity of a given directional ex- 
periment [6]. It is defined, in our case, as the minimal 
cross-section required to obtain a DM detection with a 
significance greater than the 3<r level ain 90% of the ex- 
periments. The expected sensitivity of a directional de- 
tector filled of CF 4 with a low exposure (30 kg. year) and 
a 5 keV energy threshold is displayed in fig. [TJ Such a 
directional detector should be able to reach a sensitivity 
down to a SD cross- section of 10" 5 pb for a WIMP mass 
of 20 GeV/c 2 . As one can see from the pink and cyan 
curves in fig. [T] the distance between the background 
free exclusion limit and the discovery sensitivity is lower 
at low WIMP mass than at high WIMP mass. This ef- 
fect, highlights the fact that directional detection is more 
sensitive to low WIMP mass as fewer number of events 
are required to reach a high significance detection for 
light WIMPs than for heavy WIMPs. This is of a ma- 
jor interest when considering the low WIMP mass issue 
as directional detection could bring valuable information 
to discriminate between a genuine WIMP detection and 
an unexpected background contamination of the DAMA, 
CoGent and CRESST experiments. 

For high WIMP-nucleon SD cross sections, it is also 
possible to go even further [9 . With the help of a high 
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FIG. 2: Feynman diagrams contributing to spin dependent 
elastic scattering of neutralinos on nucleons. 



dimensional mutivariate analysis, it is possible to iden- 
tify a WIMP with directional detection. It has been 
shown that dedicated analyses of simulated pseudo-data 
of a 30 kg. year CF4 directional detector would allow 
to constrain the WIMP properties, both from particle 
physics (m x ,(jp D ) and galactic halo (velocity disper- 
sions). For instance, for a benchmark model (m x = 
20 GeV/c 2 ,a^ D = 10" 4 pb), the constraints would be 
the following 

m x = 19.9±1;| GeV/c 2 (68% CL), 
logioK) = -3.97 ±0.06 (68% CL). 

Figure fl] presents the 68% contour level in the (m x ,a n ) 
plane. This is indeed a model-independent measurement 
-as the velocity dispersions are set as free parameters 
within the framework of a multivariate Gaussian veloc- 
ity distribution- of the WIMP properties, consistent with 
the input values and with a rather small dispersion. 
To assess the interest of directional detection, either 
discovery or exclusion, in the following we explore the 
MSSM and NMSSM parameter space in order to check 
if some models, excluded neither by colliders nor cosmol- 
ogy, would lie in the regions of interest. 



III. SPIN DEPENDENT ELASTIC 
SCATTERING INTERACTIONS 

A. Neutralino-nucleon spin dependent cross section 

In the non-relativistic approximation, the WIMP- 
nucleon interaction is composed of two contributions: 
the spin dependent (axial) and the spin independent one 
(scalar) [80]. Two diagrams contribute to the SD inter- 
actions of a neutralino: the Z exchange and the squark 
exchange (see Fig. [2J. The former relies on the Higgsino 
component of the LSP. When the Z contribution domi- 
nates, the amplitude for nucleon (ajy) is directly propor- 
tional to the neutralino Z coupling, 

a N = - (A? - A-d - Af ) (iV x 2 3 - ^ 2 4 ) 

where A/13, A/14 are the higgsino-d and higgsino- ?x com- 
ponents of the LSP and the coefficients describe the 
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contribution of a quark q to the spin of the nucleon. With B. Constraining the Spin Dependent cross section 



A£ = 0.842 ±0.012 A^ = -0.427 ± 0.013 
A P Q = -0.085 ±0.018 



one can show [8T] that the two amplitudes have opposite 
signs and that a p /a n — —1.14 ±0.03 independently of the 
parameters of the model. If the squark exchange is im- 
portant, the relation between amplitudes are shifted and 
can even have the same sign. Note that the relative sign 
of a p and a n is a key point for SD direct searches as it 
could lead to either constructive or destructive interfer- 
ence between the two amplitudes, depending on t he spi n 
content of the considered nucleus target (see sec. IIIB). 
The squark contribution to the amplitude reads 



a N 



t 2 w N n : 



4A 
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ru- 
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Mi) 2 (&: 



2 > 



4(m~ — ra~) 



(4) 



where Nu, N12 are respectively the bino and wino com- 
ponent of the LSP and cw = cosOw In the limit of 
an almost pure bino LSP, which is usually the case for a 
light neutralino, the amplitude simplifies to 



a N 



8M 2 c 2 



t 2 N 2 
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m~ — m~ m~ 

a FL X dpi 



(5) 



If furthermore the L and R squark masses are equal, the 
dominant contribution comes from ur and dpc respec- 
tively as the neutralino/quark/squark coupling is pro- 
portional to the quark hypercharge which is larger for 
RH quarks. The ratio of proton and neutron amplitudes 
is then a p /a n = —3.38 ±0.22. In the case where gives 
the largest contribution the ratio of amplitudes is posi- 
tive with a p /a n = 1. The squark exchange is however 
usually suppressed as compared to the Z exchange be- 
cause of the higher mass scale involved, rriq ^> Mz- It 
can be dominant when the higgsino component is very 
small leading to a suppressed coupling of the neutralino 
to the Z and to small cross sections. In fact for the squark 
contribution to be relevant for the discovery reach of di- 
rectional detection (<j p ~5x 10 -6 pb), one can estimate 
that the common squark mass has to be below roughly 
620 GeV when the LSP is a light bino (N n = 1). Such 
squark masses are excluded by LHC if the dominant de- 
cay mode of the squark is into q ± LSP [82 . Note that a 
cancellation between the squark and Z exchange can lead 
to a p ^> a n or a p <C cr n , we however stress that this oc- 
curs typically only for small cross sections if the squarks 
are near the TeV scale as will be discussed further in 
Section ED 



The Spin Dependent (SD) cross-section (at zero mo- 
mentum transfer) on a nucleus A X is given by 



r SD (A 



32 



r 2 u 2 J + 1 

U F H A -j 



(a v < S p > ± a n < S n >) 



(6) 

where Gp is the Fermi constant, ii 2 A is the WIMP-nucleus 
reduced mass, J the total nuclear spin and < S p , n > = 
< N I S Pi n I N > the neutron/proton spin content of 
the target nucleus. To be sensitive to the SD interaction, 
one may either choose a pure non-zero spin nucleus target 
( 3 He, 19 F) or rely on the isotopic fraction of the chosen 
target. This way, scalar detection experiments may also 
impose constraints on SD interaction thanks to the small 
fraction of odd-A target nucleus in the sensistive medium 
(e.g. 73 Ge for Ge-based detectors or 129 > 131 Xe for Xenon- 
based ones). 

The result of a DM experiment, either a discovery or an 
exclusion, implies a constraint on a^^X) which must 
then be converted on a constraint on a Pjfl for the sake of 
comparison between various DM experiments using dif- 
ferent target nuclei. This is however not straightforward 
as the amplitude on both proton and neutron are involved 
in ([6). 

For simplicity, the pure-nucleon coupling approximation 
(a n = or a p = 0) is often used. This method is however 
WIMP-model dependent as there is no particular reason 
for one coupling to vanish, indeed this does not occur in 
the MSSM. A model-independent method has been pro- 
posed by D. R. Tovey et al. [78 to enable comparison 
amongst SD direct searches of DM. For a given WIMP 
mass, an exclusion limit on a SD ( A X) is then translated 
into a constraint on a v and <r n , as: 



< S n > 

^ P± ^s^' 



where a l p z m is the limit on WIMP-proton cross section 
obtained in the pure proton approximation. 
The limit is expressed, for a given WIMP mass, in the 
nucleon cross- section plane (cr p ,cr n ) and divided in two 
cases: "constructive" and "destructive", whether the 
diffusion amplitude on proton and neutron add up co- 
herently or not, taking into account the relative sign 
of < S p > and < S n >. Hence, the exclusion in 
the ((jp,cr n ,m x ) space does not depend on a particular 
WIMP model. 

As outlined in [83], setting a limit on the SD interaction 
requires to neglect the SI one, which is not always justi- 
fied. In particular, in the case of Fluorine, the SD rate 
can be dominant over the SI, but this has to be checked 
in each particular WIMP model. 
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1. Proton-based versus neutron-based detectors 

The knowledge of the expectation values of the spin 
content of the proton and neutron within the nucleus 
(< Sp,n > = < N \ Sp^ n | N >) is a key issue for 
SD detection of DM. The WIMP couples mainly to the 
spin of the unpaired proton (e.g. 19 F) or to the one of the 
unpaired neutron (e.g. 3 He), leading to a contraint on ei- 
ther the proton or neutron diffusion amplitude (a ni a p ). 
One may then distinguish proton-based ( 19 F, 23 Na, 27 A1, 
35 C1, 127 I) and neutron-based ( 3 He, 73 Ge, 129 Xe) SD ex- 
periments which shall present an obvious complementar- 
ity [HUHg. 

However in practice, the spin of the target nucleus is car- 
ried both by constituent neutrons and protons and the 
comparison between SD experiments is not straightfor- 
ward. 

Detailed nuclear shell-model calculations have been de- 
velopped and the accuracy of the < S p > and < S n > 
evaluation is assessed by comparing, within the same 
shell model, the predictions on the magnetic moment and 
the energy spectra for the nuclear lowest eigenstates with 
experimental values. We refer the reader to [87 for a 
comprehensive discussion on the subject. 
Several collaborations have published exclusion limits on 
SD WIMP-nucleon cross-sections [TT HT81 123] . In the fol- 
lowing, we use (see fig. 

• for proton-based detector: COUPP-2010 [TQ, ob- 
tained with a 2-liter CF3I Bubble Chamber, with 
a 28.1 kg. day effective exposure, 

• for neutron-based detector: XENON10 [12 , ob- 
tained with 5.4 kg of fiducial liquid xenon, with 
a 136 kg. day effective exposure. 



2. The case of 19 F 

In the following we emphasize that some discrepancies 
remain on the spin content of 19 F. The value of the spin 
content of 19 F is predicted by two authors [88, 89], see 
tab. [Hj It is dominated by the proton content but the 
neutron one varies by more than one order of magnitude. 
Values from A. F. Pacheco & D. Strottman are widely 
used. However, it is noteworthy that the evaluation from 
P. C. Divari et al. is obtained by shell model calculations 
using the Wildenthal interaction [90] , which is known to 
reproduce accurately many nuclear observables. 
Nonetheless, to allow a fair comparison with existing 
published results jTTJ[l8], ^ ne vames < S p > = 0.441 and 
< S n > = —0.109 from [89] are chosen, noticing that the 
other choice would only mildly alter the result on the pro- 
ton cross-section limit but would scale up the result for 
neutron cross-section by two orders of magnitude (see. 
fig. [3]). We leave the question open, highlighting the fact 
that results on 19 F must be treated with caution owing 
to the nuclear-shell model dependence. 
However, it is worth emphasizing that, in the case of 



Model 


< s p > 


< S n > Ref. 


odd-group 


0.5 


0. 


Pacheco & Strottman 


0.441 


-0.109 [89] 


Divari et al. 


0.475 


-0.0087 [88 



TABLE II: Spin content of 19 F in various nuclear shell models. 
Henceforth, values from 89 are used. 



19 F, the spin contents having opposite sign, one expects 
a constructive interference between the proton and neu- 
tron ampli tudes, as they usually also have opposite sign 
(sec. [HIAj) 

Figure |3| presents, for a 20 GeV WIMP, the sen- 
sitivity of the considered directional detector. The 
light grey area is the region already excluded by the 
most constraining experiments (COUPP-2010 [11 and 
XENON10 [12]). In the case of XENON10, the curve 
is derived from [12] considering a pure 129 Xe detector, 
which is justified as the contribution of odd xenon 
isotopes are very close. In the case of COUPP, the curve 
is derived from [11 considering a pure 19 F detector, 
i.e. neglecting the contribution of 127 1 for such a small 
WIMP mass. The exclusion limit and the 3<r discovery 
region, for a 30 kg. year CF4 directional detector are 
presented on fig. [3] It corresponds to a rather large 
region in the parameter space, well below current limits. 



IV. TESTING SUPERSYMMETRIC 
CONFIGURATIONS 

In order to establish the reach of directional detection 
for probing supersymmetric models, we take into account 
various constraints that have an impact on the parameter 
space. 



A. Parameter space 

We consider two models, both with parameters defined 
at the electroweak scale, the MSSM and the NMSSM. 
The free parameters we take in the MSSM are the same 
as in [55]. We assume minimal flavour violation, two com- 
mon soft masses Mj and Mr for left-handed and right- 
handed sleptons, equality of the soft squark masses be- 
tween the first and second generations, Mq x 2 , while the 
mass of the third generation squarks is kept as an inde- 
pendent free parameter Mq 3 . We allow for only one non- 
zero trilinear coupling, A t . The gaugino masses Mi,M2 
and M3 are free parameters as well. In particular this al- 
lows to have Mi <C M2 , implying a light neutralino much 
below the EW scale. The ratio of the doublet Higgs VEV 
tan/5, is also a free parameter in both models. In the 
MSSM, the Higgs bilinear term, /i, and the pseudoscalar 
mass Ma are the remaining free parameters. Thus we 
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FIG. 3: Left (resp. Right) panel present the constraints in the constructive (resp. destructive) case, for a 20 GeV WIMP. The 
light grey area is the region already excluded by the most constraining experiments (COUPP-2010 [TT] and XENON10 |12|V 
For a 30 kg. year CF4 directional detector, the dark grey area is the 3a discovery region, while the solid light curve labelled 
"exclusion" is the background-free projected limit. Dashed light curve presents the same result for the alternative 19 F spin 
content values 1881. 



consider MSSM scenarios with the 11 following free pa- 
rameters 1 

Mi, M 2 , M 3 , /i, tan/3, M a ,M Il , M 1r , M~ qx 2 , M g ~ 3 , A t . 

The NMSSM is a simple extension of the MSSM with 
an additional gauge singlet superfield, S, that provides a 
solution to the naturalness problem. Indeed the param- 
eter /i = A (S) is determined by the VEV of the scalar 
singlet and is thus naturally of the EW scale [27] . In the 
NMSSM there are additional parameters related to the 
extended Higgs sector. The part of the superpotential 
involving Higgs fields reads 



W = XSH u H d 



and the soft Lagrangian is 



1 



£ soft = m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + m 2 \S\ 2 
+(\A x H u H d S + ^kA k S s + h.c). 

After using the minimization conditions of the Higgs po- 
tential, the Higgs sector, which consists of three neutral 



scalar fields, Hi, H^^ H<$ and two pseudoscalar neutral 
fields, Ai, A<i as well as a charged Higgs, H^ is described 
by six free parameters, /i, tan/3 as well as A, k, A\ and 
A K . The list of free parameters therefore contains the 
ones of the MSSM with the pseudoscalar mass, Ma, re- 
placed by 

A, ft, A\, A K , 

for a total of 14 free parameters. 

These simplified models reproduce the salient fea- 
tures of neutralino DM. Indeed, apart from the mass of 
the LSP, the most important parameters are the gaug- 
ino/higgsino content of the LSP, determined by \i and 
Mi, M2, tan/3, as well as the mass of the Higgses 
which can enhance significantly neutralino annihilations. 
Sfermion exchange, and in particular slepton exchange, 
can also play a role for light neutralinos. 

There are many similarities between the MSSM and 
the NMSSM, as will be seen in the following analysis. 
However one characteristic feature of the NMSSM is that 
the singlet fields, which mostly decouple from the SM 
fields, can be very light and yet escape LEP bounds [27] . 
Therefore it is much easier to have light neutralinos be- 
cause they can annihilate into or through the exchange 
of light singlet Higgses [91] . 



1 We do not perform a thorough exploration of the much stud- 
ied CMSSM because it is only a particular case of the MSSM. 
Furthermore in that model the LSP mass is above 50 GeV be- 
cause of the LEP limit on charginos and the relation between the 
chargino and the neutralino masses. This leaves little possibili- 
ties for parameter determination. 



B. Scanning method 

In order to thoroughly scan the parameter space we 
used a Markov Chain Monte-Carlo (MCMC) code, first 
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presented in [54 . The scanning procedure consists on 
a Metropolis-Hastings algorithm, and is based on mi- 
crOMEGAs2.4 [92H94] for the computation of all observ- 
ables. The supersymmetric spectra are calculated with 
SuSpect [95] in the MSSM and with NMSSMTools [96] in 
the NMSSM. The latter also provides collider constraints 
on the Higgs sector, on sparticles and on flavour observ- 
ables. 

Each point is generated by making a random step with 
a normal variation from the previous point in each dimen- 
sion. Then, we compute its total prior V , total likelihood 
C and total weight Q = V x C. It is kept with a prob- 
ability Mm(l, Q'/Q), where Q! is the total weight of 
the point being tested and Q is that of the source point. 
If the evaluated point is not kept, then a new point is 
generated from the last accepted point. Thus, the pa- 
rameter space is scanned via a random walk by iterating 
this procedure. 

The priors we impose are: a set of parameters has to 
lie within the boundaries of the parameter space given 



by Table III while a physical solution of the spectrum 
calculator and a neutralino LSP are required. Regarding 
likelihoods, these are displayed in Table I of [54]. We in- 
clude limits on B physics observables, on the anomalous 
magnetic moment of the muon (g — 2) M , on the Higgs 
and sparticles masses obtained from LEP and the cor- 
rections to the p parameter. In the MSSM, the limits 
on the Higgs mass were applied by making use of the 
SUSY-HIT [97] and the HiggsBounds packages [98[ [99] 
as in [55 . The HiggsBounds version used (3.1.3) include 
LEP and Tevatron results as well as first LHC results, 
more recent results from CMS presented in [100, 101 
were added a posteriori. Notice that we take the WMAP 
measurement on the DM relic density as a strict upper 
limit on the LSP relic density -obtained via the usual 
freeze out mechanism-, however, we allow the neutralino 
to have a relic density as low as 10% of the measured 
value. Indeed, the LSP could be only a fraction of the 
dark component, the rest corresponding to other dark 
particles or to a modified theory of gravity. For more 
details see [54] . 

The scans we performed in this study were aimed at 
giving a general determination of the different configura- 
tions with neutralino masses at the weak scale and below. 
However, as it was shown in [54] [55], it is difficult to find 
light (< 30 GeV) neutralinos with a random walk in both 
the MSSM and the NMSSM: the probability of falling 
in these regions that require fine-tuning is rather small. 
Indeed, in the MSSM, the neutralino LSP has to annihi- 
late via the exchange of either rather light Higgs bosons, 
scenarios that are heavily constrained by the Tevatron 
experiments as well as by CMS, or light sleptons, partic- 
ularly of staus with masses close to the LEP lower bound 
of 81.9 GeV [103]. In the case of the NMSSM, it has been 
shown in [54j 1104] that the neutralino can achieve light 
masses by annihilating via or into very light singlet-like 
Higgs bosons. 

Hence we used two different techniques to trigger the 



chains that scanned the parameter spaces. On one hand 
we let part of the chains start randomly, i.e., look ran- 
domly for a starting point with Q / 0. On the other 
hand we used the previous knowledge of fine-tuned re- 
gions explored in [54j [55] to set fixed starting points for 
the rest of the chains, in order to force the random walk 
to yield at least a few points in such regions. 

A summary of the characteristics of the runs we present 



is given in Table IV 



Parameter 


Minimum 


Maximum 


Tolerance 


Model 


Mi 


1 


1 nnn 
1UUU 


Q 
O 


both 


M 2 


100 


2000 


30 


both 


M 3 


500 


6500 


10 


both 




0.5 


1000 


0.1 


both 


tan /3 


1 


75 


0.01 


both 


M A 


1 


2000 


4 


MSSM 


A 





0.75 


0.1 


NMSSM 


K 





0.65 


0.08 


NMSSM 


A x 


-5000 


5000 


100 


NMSSM 


A. 


-5000 


5000 


100 


NMSSM 


A t 


-3000 


3000 


100 


both 


Mr 

L R 


70 


2000 


15 


both 


Mr 


70 


2000 


15 


both 




300 


2000 


14 


both 


M 9 - 3 


300 


2000 


14 


both 



TABLE III: Intervals of free parameters used for the MSSM 
and NMSSM scans (GeV units). 





Points 


Qmax 


la 


2cr 


3a 


MSSM 


1208949 


0.755 


0.25 


0.68 


0.97 


NMSSM 


2092875 


0.812 


0.30 


0.72 


0.98 



TABLE IV: Basic characteristics of the scans of the MSSM 
and the NMSSM. The la, 2a and 3 a, columns represent the 
fraction of points satisfying 0.32 x Qmax < Q, 0.05 x Qmax < 
Q and 0.003 x Qmax < Q respectively. 



V. RESULTS 

Before analyzing the predictions for SD interactions, 
we impose further constraints on the parameter space of 
the models found by the MCMC. We focus on astroparti- 
cle constraints, including XENON100 limits on SI inter- 
actions [20] as well as limits from Fermi-LAT observations 
of the photon flux from dSphs [105 2 . 



2 Here we do not show the impact of recent LHC limits on the 
MSSM Higgs sector, in particular, the LHCb and CMS combined 
results on the B s — > fifi branching ratio and the CMS negative 
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A. Dark Matter searches constraints 

Dark Matter observables are computed for each point 
kept by the MCMC. This includes the direct detection 
cross sections: SI and SD neutralino-nucleon elastic scat- 
tering processes, for both protons and neutrons, and the 
7-ray flux produced by neutralino annihilations at low 
velocities. 

Notice that the actual observables are: <J S catPo for di- 
rect detection, i.e. the elastic scattering cross section 
times the neutralino density at Earth, and cr ann pf oc for 
7-ray indirect detection, i.e. the neutralino annihilation 
cross section times the local neutralino density squared 
(in the astrophysical object). Since we allow for the neu- 
tralino to represent only a fraction of the DM component, 
we scale the neutralino density by the same fraction in 
all astrophysical systems. Thus we define 

Z = Min[ Xl 1 

\\lWMAph J 

where QwMAph 2 = 0.1097 is the la lower limit of the 
DM density as obtained by the WMAP 5-year analy- 
sis [102] . 



1. Spin Independent elastic scattering 

Spin independent interactions result from Higgs 
and/or squark exchanges. Because of the mass scales 
involved, it is generally the Higgs contribution that 
dominates, providing the LSP is some mixture of hig- 
gsino/gaugino. In the MSSM, both the light and heavy 
Higgses can contribute, in particular the heavy Higgs 
contribution is enhanced at large values of tan (3. In the 
NMSSM, there is in addition a contribution from the sin- 
glet Higgs. When the lightest scalar is below 10 GeV 
one can have an enhancement of the cross section even 
if the LSP is weakly coupled to the scalar Higgs j!04j . 
The SI interaction rates depend on the quark content of 
the nucleons. Indeed, the choice of quark coefficients 
in nucleons could vary the estimation by an order of 
magnitude. In micrOMEGAs the quark content of nu- 
cleons is parametrized by the a^N and <Jo terms [92] . 
Recent lattice QCD results point towards small strange 
quark contributions, hence to a n N — 0"o 1107] . We take 
a^N = 45 MeV and a = 40 MeV. 

The predictions for SI interactions in the MSSM and 
the NMSSM along with XENON100 [20^ and CDMS- 
II [23] limits are displayed in Fig. [4] Color tagging in 
relation with detectability in directional detection will 
be discussed in Sec. IVCl The XENON100 limit attain 



search for H — > tt constraining the tan (3 vs. Ma plane. For 
an account of such an analysis over the same MSSM data set 
presented here, see 106 . 




FIG. 4: Spin independent cross section versus the neutralino 
mass. Top: MSSM. Bottom: NMSSM. In pink the points in 
the discovery region of directional detectors and in cyan in 
the exclusion region. The CDMS-II [23] and XENON100 [20] 
limits are also displayed. 



some of the configurations found in both models, thus 
constraining the parameter spaces. Nevertheless, most 
scenarios lie below the XENON100 sensitivity, and many 
are more than an order of magnitude away from it. This 
highlights the need for complementarity in DM search 
strategies. 

While in general the point distribution is rather sim- 
ilar in both models, there are a few differences between 
the two panels. Those differences can best be understood 
by investigating the effet of the relic density constraint 
which depends on the model. In particular, neutrali- 
nos below 30 GeV in the NMSSM usually imply a light 
Higgs scalar or pseudoscalar to ensure the annihilation 
rate is large enough to obtain relic densities below the 
WMAP value. If the light Higgs is a scalar, then the 
SI cross section gets enhanced, as shown in |104[ 1108] , 
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whereas a light pseudoscalar does not contribute to the SI 
cross section. This explains the larger range of predicted 
cross sections observed in the NMSSM with respect to 
the MSSM. In the latter, as discussed in [55], when 
the relic density is achieved through as-light-as-possible 
scalar Higgs exchanges, the lighter the neutralino, the 
wider the mass difference between neutralinos and Hig- 
gses, thus the larger the coupling between these must 
be. Hence the SI cross section tends to increase towards 
lighter neutralino masses, which corresponds to the top- 
left arm of the cloud in Figs. [4js top panel. However, 
when the neutralino annihilations are driven by slepton 
exchanges, the SI cross section does not evolve with the 
neutralino mass, and is not necessarily enhanced by the 
need of large neutralino-Higgs couplings. That situation 
corresponds to the other light neutralino arm, a more dif- 
fuse cloud at smaller interaction rates, in the same panel. 

Above 30 GeV, neutralinos can achieve the correct 
relic density by resonant Z exchange. That region (for 
m % o ~ 45 GeV) is well represented in both models. In 
both cases many points show small SI rates. Since the 
Z resonance enhances the annihilation rate, the ZxiXi 
coupling should be suppressed, in order to have enough 
neutralinos left after the thermal freeze-out. Hence the 
Higgs couplings, which depend on the higgsino compo- 
nent of the LSP as the Z coupling, are also suppressed, 
which in turn diminishes the SI cross section. Above 
the Z resonance, and around 60 GeV, the SM-like Higgs 
exchange on resonance dominates the relic density. For 
the same reason, when the masses are very fine-tuned, 
m % o & rrih/2, the required annihilation rate forces the 
Higgs-neutralino coupling (responsible for the interac- 
tion) to be fairly small, therefore allowing very weak SI 
interactions. Note that the 45-85 GeV neutralino mass 
range with large cross-sections is more populated in the 
MSSM than the NMSSM. This could be understood in 
terms of the couplings of the SM-like Higgs to the LSP 
which can be suppressed because of the singlet compo- 
nent in the NMSSM. 

Approaching 100 GeV, neutralinos could have larger 
higgsino components, since we can have fi > 100GeV ~ 
m % o. Furthermore there are many possibilities for neu- 
tralinos to achieve the relic density: by Z or Higgs ex- 
changes (large higgsino components) or via sfermion and 
gaugino exchanges. Also, the heavier the neutralino, the 
narrower the neutralino-squark mass difference could be. 
This enchances the squark contribution and allows for 
some unusually large cross sections. The predictions for 
the SI cross section span five orders of magnitute. In 
general, the heavier the neutralino, the smaller the cross 
section and the XENON 100 limits eventually do not con- 
strain many configurations with m % o < 200 GeV. 



2. Gamma rays 

Neutralinos in galactic objects, such as the Milky Way 
and its dwarf spheroidal companions (hereafter referred 



to as dSph), have a probability of encountering and anni- 
hilating into SM fermions. After the subsequent decays 
and hadronization, these events produce 7-rays. More 
marginally, 7-rays can be produced from internal lines 
in the annihilation process or directly by pairs, the lat- 
ter which occur through a loop-induced process is not 
included in our flux computation. Among the indirect 
signatures, 7 rays can be easily computed, as compared 
to charged cosmic rays, they do not suffer from energy 
losses and their propagation is not subject to uncertain- 
ties on the interstellar medium. Also, the Fermi-LAT, 
HESS and MAGIC experiments are successfully explor- 
ing the fluxes at Earth, thus providing means to constrain 
DM models. 

The expected flux of 7 photons of energy E for a given 
angular spread iJj : (j) 1 {E,ijj), depends upon the square 
of the neutralino density, on the annihilation cross sec- 
tion (<jv) times the squared local abundance of neutrali- 



nos p and on the particle spectrum produced 



dN(E) 



"X? *" *" *" dE - 

The computation of the flux implies the integral over the 
line of sight (l.o.s.) I within the extension of the desired 
area, which is usually matched to that of an observatory. 
Hence 



-J 
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di^) P io(im x 



1 (<jv) dN(E) 

2 47rm 2 dE 



We have explicitly split the so called astrophysics term 
J (ip) -namely the l.o.s. integral containing all the terms 
which carry spatial dependence- and particle physics 
term <p pp (E) -which is the energy dependent term. The 
former can only be determined with the knowledge of the 
DM distribution in the observed object. We focus here 
on dSphs observed by Fermi-LAT |105j . for which they 
provide an estimation of J (ijj) (see their Table 4). We 

compute J (j) PP (E) over the 100 MeV, m % o interval. 

Here, we have scaled down the flux by a £ 2 factor in or- 
der to take into account only the neutralino contribution 
to the DM density distribution. 

The annihilation cross section at galactic velocities is 
related to the interaction rate at freeze out. However, 
the velocity in galaxies is much lower than at freeze out, 
hence large variations can take place when the annihi- 
lation occurs via resonances in the early universe [109J. 
In general annihilations that proceed through a Higgs 
and/or Z resonance in the early universe are diminished 
at v ~ 10 _3 c, furthermore sfermion exchange into light 
fermions get suppressed. However it is also possible to 
have an enhancement of the neutralino annihilation cross 
section in galaxies when it is dominated by the exchange 
of a Higgs with a narrow width, see [104 for a discus- 
sion on this subject in the particular case of the light 
neutralinos in the NMSSM. 

Fig. [5] shows the 7-ray integrated flux expected from 
the Draco dSph as a function of the neutralino mass, 
along with the Fermi-LAT limits [105 . Color tagging in 
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FIG. 5: Flux of 7-rays expected from neutralino annihilations 
from the Draco dwarf spheroidal galaxy versus the neutralino 
mass. Top: MSSM. Bottom: NMSSM. In pink the points 
in the discovery region of directional detectors and in cyan 
in the exclusion region. The Fermi-LAT limits [105] are also 
displayed. 



FIG. 6: Proton-neutralino spin dependent elastic scattering 
cross section versus the neutralino mass with the exclusion 
and discovery projections for a nominal directional detector. 
Top: MSSM. Bottom: NMSSM. In blue safe points and in 
yellow points excluded by either XENON100 or Fermi-LAT. 



relation with detectability in directional detection will 
be discussed in Sec. |V C| Fewer points are challenged 
by Fermi-LAT than by XENON100. Only the lighter 
configurations can be constrained by their 7-ray yield, 
while heavier neutralinos seem to be out of the reach 
of the Fermi-LAT detector. As for SI interactions, re- 
sults for the MSSM and the NMSSM differ mainly for 
the lighter neutralino configurations. Here, the scalar 
Higgs exchanges suffer from a low-velocity suppression, 
while the pseudoscalar exchanges can be resonantly en- 
hanced, this is possible only in the NMSSM. Thus, again, 
the NMSSM predicts a broader range of possible 7-ray 
yields when neutralinos are light. An order of magnitude 
enhancement of the Fermi-LAT sensitivity would probe 
all the configurations in the MSSM up to 30 GeV neu- 



tralino masses, and many scenarios up to 200 GeV in 
both models. 



B. Predictions for Spin Dependent interactions 

The predictions for the £(Tp D observable, which is the 
target of the future directional detectors, can be seen in 
Fig. [6] We also show the projected curves for the dis- 
covery limit (in pink), defined as a significance greater 
than the 3<r level at 90% Confidence Level (C.L.) in 30 
kg. year CF4 directional detector, and the exclusion limit 
(in cyan), in the background- free case for the same detec- 
tor with the same exposure, see sec. [TTJ These curves de- 
fine three regions: the discovery region above the discov- 
ery curve, the exclusion region between the two curves, 
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FIG. 7: Spin dependent elastic scattering cross sections cor- 
relations: proton-neutralino versus neutron-neutralino inter- 
actions. Contrarily to fig. [3] the result is presented for all 
neutralino masses. Top: MSSM. Bottom: NMSSM. In pink 
the points in the discovery region of directional detectors and 
in cyan in the exclusion region. 

and the out-of-reach region below the exclusion limit. We 
have tagged points failing to overcome the XENON 100 
and/or the Fermi-LAT constraints in yellow. MSSM and 
NMSSM configurations range over several orders of mag- 
nitude, the maximum reaches 10 -39 cm 2 and a large frac- 
tion of the points lie above the potential exclusion limit 
of future directional detectors. One can readily see that 
some configurations that are not yet constrained by SI 
detectors nor by indirect signals (blue points) lie in the 
exclusion or even the discovery region. 

In particular most neutralinos with a mass m % o < 40 
GeV lead to large cross sections. This is understand- 
able since in this case the same diagram (Z exchange) is 
responsible for both annihilation and scattering on nu- 
cleon. Enhancing the annihilation cross section in order 
to get an acceptable relic density implies enhancing the 
SD cross sections. However, when the neutralino mass 
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FIG. 8: Frequency distribution of ap/an. Top: MSSM. Bot- 
tom: NMSSM. In pink the points in the discovery region of 
directional detectors and in cyan in the exclusion region. 



approaches mz/2, the Z-neutralino coupling is reduced 
so the relic density is not too small, hence yielding very 
small SD interaction rates. 

For larger neutralino masses, the possible cross section 
values range similarly to the SI case. However, here there 
is a clear upper limit for the SD, Z-exchange preferred 
configurations. For m % o > 100 GeV the SD cross section 
scales as l/m 2 . Also, the largest coupling between 

A i 

the Z and the neutralino is Max(N 2 3 — 7V 2 4 ) = 0.5, 
which can be met for higgsino dominated neutralinos. 
Those configurations having the maximum coupling 
describe the upper limit of the cloud above 70 GeV, 
drawing the l/m 2 curve. Points falling above that limit 
imply dominant squark exchanges, and points below 
have either smaller couplings or destructive interference 
between the Z and squark exchange diagrams. 
It is worth emphasizing that in the MSSM, and to a 
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lesser extent in the NMSSM, a large fraction of the 
supersymmetric configurations below m % o < 200 GeV 
lies in the discovery region, meaning that such models 
could be discovered, with a significance greater than 3a 
(90% CL), with a 30 kg. year CF 4 directional detector. 
Exclusions may be reached up to ~ 800 GeV. 

As outlined in sec. |III[ a model independent con- 
straint on DM requires to present the results in the 
((Tp D ,(jS D i m x) parameter space 3 . In order to have a 
complete view of the theoretical predictions, we present 
on Figs. [7] the cr SD on proton versus neutron, for all neu- 
tralino masses. The discovery (resp. exclusion) limit 
depends on the mass, ranging from ~2x 10 -41 cm 2 (~ 
7 x IO- 42 cm 2 ) for a - 10 GeV WIMP to - 2 x lO" 40 cm 2 
(~ 2x 10- 41 cm 2 )) for a - 500 GeV WIMP. Color tagging 
refers to the detectability in directional detection: pink 
discovery region, cyan exclusion region and blue out-of- 
reach. As expected, in models domin ated by Z exchange, 
we found cr p /a n = 1.3 (see Sec. Ill A). Most of the points 



in the discovery region (pink points) satisfy this condi- 
tion. Exceptions are points where Nf 3 — Nf 4 is small and 
the squark are light (more specifically below 300 GeV). 
The latter could easily be excluded by the LHC, unless 



they are degenerated with neutralinos, see Sec. [VI] for a 
more complete discussion on that respect. 

When the squark exchange comes into play, inter- 
ference between the squark and Z contributions can 
reduce either the neutron or proton cross sections, and 
the correlation is lost. Indeed, when the Z exchange 
is dominant but the squark exchange becomes non 
negligible, the proton cross section is lessened rather 
than the neutron cross section (points deviating from 
the correlation line towards the bottom). When the Z 
and squark exchanges are similar, then the destructive 
interference can completely erase one or the other of 
cross section, but not both, depending on the nature 
of the lightest squark. These cases correspond to 
the broad vertical (proton cross section suppression) 
and horizontal (neutron cross section suppression) 
distributions towards smaller cross sections in the top 
panel of Fig. [7] Finally, when the squark exchange 
dominates, the proton cross section being fa vored, the 
correlation reaches <J p /<J n = 11.4 (see Sec. Ill A). In 
the NMSSM the correlation <J p /<J n = 1.3 extends to 
lower values of the cross-section. In such a case the 
neutralino is dominantly higgsino with a significant 
singlino component, furthermore N13 « N14 so that 
the Z contribution is suppressed, hence the low cross 
section. The annihilation of neutralinos in the early 
Universe is driven by the singlino component and the 
Higgs sector. Such scenarios are more difficult to obtain 
in the MSSM, furthermore the lower bound on the relic 



density constrains the case of a dominantly higgsino LSP. 

Fig. [8] presents the frequency distribution of a p /a n for 
the points of Fig. [7] with the same color code. As outlined 
in sec. |IIIB the relative sign of a p and a n is a key issue 
for SD direct searches. Indeed the interference can be 
either constructive or destructive, depending of the sign 
of the spin contents of the target nucleus (see eq. [6J. For 
both the MSSM and the NMSSM we observe that most 
models give a p /a n ~ —1.13, corresponding to a dominant 
Z exchange. A small number of models give a p /a n ~ 
—3.5, corresponding to a squark (qn) contribution in the 
pure bino limit. However, amplitudes can also be of the 
same sign. As we shall see in section VI, this is not the 



case when applying recent ATLAS limits on squarks and 
gluinos. 



C. 



The complementarity between directional 
detectors and other techniques 



3 Note however that it requires to set the SI coupling to zero, which 
is not model independent stricto sensu. 



We have already shown in Sec. VB that directional de- 
tectors will scan an important fraction of the supersym- 
metric parameter space. We have also shown the power 
of XENON100 (in direct SI detection) and Fermi-LAT 
(in indirect detection experiments) to put constraints 
on neutralino DM configurations in Sec. VA For the 
NMSSM case, it has been discussed in |104j that SI di- 
rect detection and indirect detection of 7-rays constrain 
different light neutralino scenarios. In particular the light 
scalar Higgs case is reachable by the former, while light 
pseudoscalar configurations are probed by the latter. 

Let us discuss the impact of SD direct detection 
searches and their position in this sense. Figs. [4] and [5] 
show the yields in SI cross sections and in 7-rays from 
the Draco dSph respectively. In these planes we have 
tagged in pink the scenarios that fall in the discovery re- 
gion of future directional detectors, in cyan scenarios be- 
tween the exclusion and the discovery curves, and in blue 
those points which are not expected to produce enough 
signal. In both cases one can see that there are many 
pink points below the exclusion limits of XENON100 (in 
Figs. |4]) and below those of Fermi-LAT (in Figs. [5]). Most 
of these pink points correspond to m % o < 200 GeV. We 
can observe that the development of a SD~ directional de- 
tector will help closing the parameter space that seems 
to escape other techniques. Unfortunately, we also ob- 
serve many blue points in both figures, which lie orders 
of magnitude below the reach of all DM search strategies. 

We summarize the complementarity of detection tech- 
niques in Figs. [9j in which we show the correlation be- 
tween SI and indirect detection techniques, using the 
color tagging for the directional detectors potential reach. 
Those points failing to satisfy the XENON100 limit 
(which cannot be represented in this plane) are not 
drawn. Notice that the pink cloud, set to be upfront, 
covers a large region of this correlation plane, spanning 
over more than 2 (7) orders of magnitude in the £,cr SI axis, 
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FIG. 9: Flux of 7-rays expected from neutralino annihilations 
from the Draco dwarf spheroidal galaxy versus spin indepen- 
dent cross section. Top: MSSM. Bottom: NMSSM. In pink 
the points in the discovery region of directional detectors and 
in cyan in the exclusion region. The Fermi-LAT limits 105 
are also displayed. Points excluded by XENON 100 are not 
shown. 



and over more than 5 (7) orders of magnitude in the <p 1 
axis for the MSSM (NMSSM). This, of course, includes 
many configurations that are not only far away from the 
current reach of SI and indirect detectors, but that would 
also escape their forthcoming upgrades. Supersymmetric 
models with neutralino DM candidates predict possible 
configurations in the range of detection for all the exper- 
imental techniques exposed here. However, it is obvious 
that none is able by itself to probe every configuration we 
found. Indeed, the reach of any of these techniques alone 
is frustratingly insufficient. In this sense, the develop- 
ment of SD detection techniques is useful for a thorough 
scan of neutralino DM supersymmetric configurations. In 
particular, the superposition of the blue, cyan and pink 
clouds in Figs. [9] leads us to conclude that the capabil- 
ity of SD directional detection, either to discover or to 
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FIG. 10: Correlation plot between the fi mass term and the 
bino mass Mi. Top: MSSM. Bottom: NMSSM. In pink the 
points in the discovery region of directional detectors and in 
cyan in the exclusion region. 



exclude a DM model, is not correlated with the sensitiv- 
ity of the SI and indirect searches. It does in particular 
highlight the complementarity of the various DM search 
strategies. 



D. Constraining MSSM and NMSSM parameter 
spaces with directional detectors 

Although the MSSM and NMSSM parameter spaces 
are multidimensional and can produce many different 
phenomenologies, the reach of directional detectors is de- 
termined solely by the DM particle mass and the SD cross 
section. In particular, for masses below 200 GeV, a large 
fraction of the supersymmetric configurations fall in the 
discovery region. For neutralinos to satisfy m % o < 200 
GeV, the mass term of at least one of its components 
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FIG. 11: Free parameter frequency distribution normalized to 
Qmax in the MSSM run. In pink the points in the discovery 
region of directional detectors and in cyan in the exclusion 
region. 





has to satisfy the same condition. This can be either the 
bino component determined by Mi, the higgsino com- 
ponent given by /i, the wino component given by 
or the singlino component in the NMSSM (which de- 
pends on a combination of A, k and n). Note however 
that when the neutralino LSP is wino-dominated, which 
implies < Mi,/i, the relic density tends to be very 
small, those configurations are therefore excluded by the 
condition & x oh 2 > 10%Q\y m Aph 2 . The neutralino LSP 
main component is therefore determined by the lighter of 
Mi,/i. The strength of the SD interactions, strongly de- 
pend on the \i value, determining the neutralino mixing 
elements 7Vi 3 and 7Vi 4 , thus the ZxiX? coupling. Since 
the lighter the /i the larger | 7Vi3 | and | 7Vi4 |, the 
lighter the neutralino (although depending on Mi, and 
the singlino component in the NMSSM), we expect small 



FIG. 12: Free parameter frequency distribution normalized to 
Qmax in the NMSSM run. In pink the points in the discovery 
region of directional detectors and in cyan in the exclusion 
region. 



\i values to be generally in the reach of directional detec- 
tors. 



In Fig. [10] we present the correlation between [i and 
Mi, using the same color tag for points in or out the 
discovery and exclusion regions. The general trend just 
discussed is indeed met: a discovery is mostly linked with 
\i < 250-300 GeV. Roughly speaking, the \i = Mi line 
separates the bino and higgsino dominated neutralino 
compositions. 

For a bino-like neutralino to have large enough SD in- 



16 



teractions (pink or cyan points), either the higgsino com- 
ponent is not negligible (i.e. points above the /i = Mi line 
but not far away from it), either the squarks have similar 
masses to those of neutralinos (which is the case of the 
few pink and cyan points above and far away from the 
H = Mi line). The former case is very well represented 
by the vertical clouds at Mi ~ 50 GeV, corresponding 
to annihilation via the Z resonance. In such a line, \i 
below 300 GeV allows discovery (pink points), and up 
to ~ 450 GeV the configurations could be challenged by 
directional detectors (cyan points). Larger values of \i 
would be out of reach (blue points). 

The effect of the neutralino mass is imprinted in the 
fact that very few pink points can be found at large values 
of Mi and \i. This is easily understood: the heavier the 
neutralino, the smaller the cross section and the smaller 
the projected sensitivity, the less configurations fall in 
the discovery region. 

The correlation between small values of \i and the dis- 
covery potential is also observed in Figs. 11 and 12 In 



these figures we show the frequency distribution of all 
free parameters in the MSSM and NMSSM respectively 
for all points (blue), as well as for pink and cyan points 
only. In the fourth panel of both figures one can see that 
the pink distribution of /i is strongly peaked towards the 
lighter values, almost reaching the all point curve in the 
MSSM. Indeed, few configurations with \i < 150 GeV 
would actually escape a discovery with a large directional 
detector (30 kg. year). In the NMSSM such configurations 
could also have a large singlino component in which case 
they would remain out of reach of directional detectors. 
Regarding M x (first panel), the pink and cyan distribu- 
tions also reflect the neutralino mass range for detection 
and exclusion. In particular, apart from some configu- 
rations with Mi around Mz/2 (Z resonant neutralinos), 
most configurations with Mi < 450 GeV would be chal- 
lenged by future directional detectors. Conversely, for 

Mi > 450 GeV, most supersymmetric configuration 
would not be probed by directional detectors, even if ex- 
clusion remains possible. 

As for the other parameters, the shapes of distribu- 
tions are to a large extent related to the LSP mass. The 
exclusion region does not contain any points with a LSP 
lighter than 35 GeV and extends to much larger masses 
for the LSP than the discovery region. Thus the cyan 
distributions for the masses of supersymmetric particles 
will be shifted towards higher values than the pink dis- 
tributions. In particular the distributions for the soft 
terms /i, M2 that drive the masses of the neutralinos and 
charginos are peaked at higher values for the exclusion 
region than for the discovery region. For the same rea- 



son, the pink distribution for Mr 



and M r 



are shifted 



towards smaller values than the cyan and blue distribu- 
tions. Regarding tan /3, we also expect to see that smaller 
values are preferred for the points in the discovery re- 
gion. Indeed, in the MSSM, the light neutralino LSP 
is often associated with a light M^, collider constraints 
on the Higgs sector then selects small values of tan/3. 



This trend is also met in the NMSSM, however, since the 
light neutralino is not necessarily associated with a light 
pseudoscalar mass, there is a much larger fraction of con- 
figurations with small values of tan /3 that escape SD de- 
tectability. Finally we remark that the masses of squarks 
and gluinos can extend well above 1 TeV, this sector is 
mostly constrained by flavour observables. Therefore a 
large fraction of the configurations will escape detection 
at the LHC as will be discussed next. 



VI. IMPACT OF LHC RESULTS ON THE 
PARAMETER SPACE 
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FIG. 13: Top Panel: Spin dependent elastic scattering 
cross sections correlations: proton-neutralino versus neutron- 
neutralino interactions in the MSSM. In pink the points in 
the discovery region of directional detectors and in cyan in 
the exclusion region. Here, in contrast with Figs. [7| we have 
removed those points falling above the ATLAS limit on the 
Mq vs. Mg plane. Bottom panel: frequency distribution of 
a p /a n for the same points with the same color code. 



17 



The ATLAS collaboration has published an analysis 
and established a limit on the squark mass versus gluino 
mass plane [82]. These limits were derived for mass- 
less neutralino LSP in the MSSM and assuming that 
the squarks decay exclusively into a quark and the LSP. 
This is not the case in the scans we have performed in 
the MSSM, in particular additional decay modes of the 
squarks can weaken the ATLAS limits. However, com- 
paring these limits to our data sets is a good way to 
check the possible influence of the forthcoming results of 
the LHC in terms of parameter space cutting and its im- 
plications for SD cross sections. The modification of the 
squark decay modes and the presence of light Higgs states 
mean that these limits do not apply to the NMSSM, thus 
we restrict this analysis to the MSSM 4 . 




m x0 [GeV] 



We have applied the limit on the Mq vs. Mg plane, 
where Mq stands for the lightest first or second genera- 
tion squark mass. Points falling below the curve provided 
by ATLAS are not represented in Fig. 13 (top panel), 
where we show the correlation between neutralino SD in- 
teractions with protons and neutrons. The difference is 
quite striking with respect to the top panel in Fig. [7| 
where all points are drawn. What we observe is that, as 
the LHC probes the lightest squarks (in the hypothesis 
these are not observed), only the points in the main cor- 
relation line are left. Thus, the possibilities to achieve a 
large SD interaction would be restricted to a large hig- 
gsino fraction. 



Bottom panel of Fig. [13] presents the frequency distribu- 
tion of a p /a n for the same points with the same color 
code. As outlined in sec. |III B the relative sign of a p 
and a n is a key issue for SD direct searches. It is worth 
noticing that taking into account recent ATLAS limits 
on squarks and gluinos restricts the allowed models to 
the ones with negative values of a p /a n . More precisely, 
only the Z exchange is allowed, and the ratio o f the two 
amplitudes is given by a p /a n ~ —1.14 (see sec. Ill A). It 
implies that a nucleus target with spin contents of oppo- 
site sign (e.g. 19 F, 3 He or 133 Cs) will present a construc- 
tive interference while nucleus target with spin contents 
of same sign (e.g. 129 Xe, 131 Xe or 73 Ge) will have a de- 
structive interference thus reducing the cross section on 
nucleus and hence the event rate. A dedicated study is 
needed on this issue, but we emphasize that it may have 
consequences in the choice of target for upcoming SD 
experiments. Note that this statement is rather indepen- 
dent of the details of the SUSY model, as long as LHC 
exclusions on squarks apply. 
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FIG. 14: Proton-neutralino spin dependent elastic scattering 
cross section versus the neutralino mass in the discovery sce- 
nario runs. Top: MSSM. Bottom: NMSSM. We display the 
la and 3a contours, while we used the former as the con- 
straint for the random walk. Also displayed are the exclusion 
and discovery projections for a nominal directional detector. 
We display only safe points regarding XENON 100, Fermi- 
LAT and CMS (MSSM only) limits on SI elastic scattering 
interactions, 7-rays from the Draco dSph, and Higgs interac- 
tions. 



VII. A DISCOVERY SCENARIO 



4 An analysis of the ATLAS constraints on the NMSSM parameter 
space is ongoing 



As outlined in sec. |IIB 2[ a dedicated analysis of data 
of a 30 kg. year CF4 directional detector could also allow 
us to constrain the WIMP properties, both from parti- 
cle physics (m Xl (j^ D ) and galactic halo (velocity disper- 
sions) [9 . This would be one step beyond current DM 
search strategy capabilities. Of course, this requires a 
rather large SD cross section associated with a low neu- 
tralino mass. The outcome would be a constraint on the 
mass and the cross section, i.e. within a contour defined 
by a confidence level. 
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FIG. 15: Free parameter frequency distribution normalized 
to Qmax in the MSSM discovery scenario run. 



In light of such a tool, we perform the following exer- 
cise [9]: we assume the existence of a DM particle of 20 
GeV mass with 10 -4 pb interaction rate with protons, 
leading to w 80 WIMP events in a 30 kg. year CF4 di- 
rectional detector, to generate simulated data. Then, we 
infer the contour from the data analysis procedure after 
detection by a canonical directional detector. We exploit 
this result by including a new prior in the MCMC de- 
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FIG. 16: Flux of 7-rays expected from neutralino annihila- 
tions from the Draco dwarf spheroidal galaxy versus spin in- 
dependent cross section. Top: MSSM. Bottom: NMSSM. The 
Fermi-LAT limits [105 are also displayed. Only the points 
overcoming all constraints (from XENON100, Fermi-LAT and 
CMS) in the discovery scenario are displayed. 



scribed in IV B we impose the neutralino mass and the 
SD cross section to lie in the la contour obtained by 
the discovery. Hence, we are able to scan the possible 
scenarios fitting the observation. 



In Fig. 14 we show the distribution of points within 
the contour in the SD vs. neutralino mass plane. Notice 
that the transition in the contour prior was sharp, so no 
point is kept outside those boundaries. In the NMSSM 
(bottom panel), we have much less statistics than in the 
MSSM (top panel), however, populating more this dis- 
covery region is only a matter of time. Nevertheless, 
the contour is much more homogeneously filled in the 
NMSSM case, showing that there are more possibilities 
to find good configurations, in contrast with the MSSM 
were the heavy end of the contour is much more preferred 
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than the actual value of the WIMP in the simulation in- 
put. Hence, the NMSSM neutralino DM appears a much 
more plausible explanation if such a detection was made. 
In these figures we have kept only the points that sat- 
isfied all constraints including the CMS limits on tan/3 
vs. M A plane (MSSM only) as well as Fermi-LAT and 
XENON100 limits. 



An observation would allow to refine the expectations 
on the soft SUSY parameters, as can be seen by compar- 



ing the parameter distribution in Fig. 11 with the one cor- 



responding to points in the contour in Fig. 15 in the case 
of the MSSM. Most distributions are determined by the 
condition on the neutralino LSP mass. This is clearly the 
case for M\ which drives the mass of the light bino, but 
also applies to other parameters [55]. Indeed for a light 
LSP to be consistent with the relic density constraint 
requires either a light slepton, hence the very peaked dis- 
tribution in for Mj r or a light Higgs boson, hence the 
peak in the Ma distribution. Higgs and flavor physics 
constraints then imply that intermediate values of tan /3 
are disfavored. Furthermore the squark contribution that 
is needed to cancel the Higgs contribution in B(b — >• sj) 
explains the peak at low values of the third generation 
squark masses Mq 3 . The main impact of imposing a spe- 
cific range for a SD is reflected in the distribution for /i, 
confined to be in a narrow range. As mentioned previ- 
ously, the strength of the neutralino coupling to the Z is 
the most relevant parameter in computing cr SD . 



Figures 16 represent the safe points in the 7-rays from 
Draco vs. SI interaction plane. There is a concentra- 



tion of points at 10 46 

-13 ^ vv ,3 -i 



10 



-45 



cm cross sections and 



10 -15 cm 5 s _i fluxes. This constitutes the prediction 
for other detection techniques in case an observation is 
made. Due to the interplay of parameters, the spread of 
these predictions is rather large in the NMSSM: 9 orders 
of magnitude in the SI axis and 6 orders of magnitude 
in <p 1 . An order of magnitude spread is predicted in the 
MSSM. This is an extremely different behavior of the 
two models. Such a discovery would predict a MSSM 
neutralino to be found shortly by direct detection and 
indirect detection experiments, thus provinding a cross 
check of the discovery claim. On the contrary, a NMSSM 
neutralino could easily escape any other detection. 



In any case, keep in mind that such a discovery would 
still mean that \i < 200 GeV, and 10 GeV < m % o < 
30 GeV. With these two characteristics we can predict 

the presence of a chargino with m + < 200 GeV. A light 

x 1 

slepton would be favored in the MSSM case, but its ab- 
sence leaves only the NMSSM as a possibility. With the 
help of collider physics, it could be possible to rule out 
the possibility of a MSSM neutralino DM, or even of any 
supersymmetric neutralino if no charged particle is ob- 
served in the vicinity of the weak scale. 



VIII. CONCLUSION 

We have shown that supersymmetric models with neu- 
tralino DM predict signals in the range of detection for 
SI and SD direct detectors as well as for indirect detec- 
tors. However only a fraction of the parameter space of 
either the MSSM or the NMSSM can be probe by each 
type of experiments alone. The development of SD di- 
rectional detection techniques thus offer the possibility 
to probe more thoroughly the parameter space of super- 
symmetric models. This is the case even if the sensitivity 
of indirect detectors as well as SI detectors is increased 
by one order (or more) of magnitude, thus emphasizing 
the complementarity between different techniques. 

With the planned MIMAC detector, neutralinos up to 
200 GeV could be discovered and up to 600 GeV could 
be excluded. The light neutralinos that are best probed 
by directional detectors are often accompanied by not so 
heavy charginos/neutralinos and even sleptons or by a 
light pseudoscalar. The search for these weakly inter- 
acting particles at the LHC will therefore impact in the 
future the potential of directional detection to probe su- 
persymmetric models. We expect the LHC to consider- 
ably expand the constraints on supersymmetric scenarios 
in the next year, even if no signal of physics beyond the 
standard model is found. We have shown that if squarks 
of the first two generations are excluded up to a mass 
of nearly 750 GeV, notwithstanding that they could es- 
cape detection due to small mass splittings with some 
other supersymmetric particles, the squarks play little 
role in direct detection and the SD cross section on pro- 
tons is completely correlated with that on neutrons. The 
LHC is also probing supersymmetric models with Higgs 
searches, in particular the negative search results on the 
heavy Higgs doublet of the MSSM constrain the super- 
symmetric models with a light neutralino. Furthermore a 
confirmation of a Higgs signal at 125 GeV as announced 
by ATLAS would narrow down the number of allowed 
supersymmetric configurations. Note however that the 
light Higgs is not directly linked with SD direct detec- 
tion. 
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